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Wc characterize all fixed equilibrium point singularity distributions in the plane of 
logarithmic type, allowing for real, imaginary, or complex singularity strengths T. 
The dynamical system follows from the assumption that each of the N singulari- 
ties moves according to the flowfield generated by all the others at that point. For 
strength vector T e M^, the dynamical system is the classical point vortex system 
obtained from a singular discrete representation of the vorticity field from incom- 
pressible fluid flow. When F G Jm, it corresponds to a system of sources and sinks, 
whereas when F e the system consists of spiral sources and sinks discussed in 
Kochin et. al. (1964). Wc formulate the cquilibrimn problem as one in linear alge- 
bra, ^F = 0, A e C^^^, F e C^, where A is a iV x complex skew-symmetric 
configuration matrix which encodes the geometry of the system of interacting sin- 
gularities. For an equilibrium to exist, A must have a kernel. F must then be an 
element of the nuUspace of A. We prove that when N is odd, A always has a kernel, 
hence there is a choice of F for which the system is a fixed equilibrium. When N 
is even, there may or may not be a non-trivial nullspace of A, depending on the 
relative position of the points in the plane. We describe a method for classifying the 
equilibria in terms of the distribution of the non-zero eigenvalues (singular values) 
of A, or equivalently, the non-zero eigenvalues of the associated covariance matrix 
A^ A, from which one can calculate the Shannon entropy of the configiiration. 

Keywords: Singularity dynamics; Fixed equilibria; Singular values; Point 
vortex equilibria; Shannon entropy 



1. Introduction 

Consider the vector field at ^ = governed by the complex dynamical system: 

z* = ^^, z{t)eC, Fee, teR>0, (i.i) 

ZTTl Z 

where z* denotes the complex conjugate of z{t). Letting z{t) = r{t)ex.p{iO{t)), 
T = Tr + iTi, gives: 
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from which it is easy to see that: 

rit) = 



27r 



9{t) 



^'■^lnii^)t + r^{0)] ifr, ^0 



(1.4) 



(1.5) 



27rr2(0) 



6(0) 



if r,, = 0. 



When Tr ^ 0, Ti = 0, the field is that of a classical point- vortex (figure [ij a), (b)); 
when Tr = 0, Ti ^ it is a source (r, > 0) or sink (r, < 0) (figure [l];c),(d)), while 
when 7^ 0, F,; ^ 0, it is a spiral-source or sink ( (figure [ije)-(h)). 







(a) r,. > 0, Ti = (b) < 0, Ti = (c) = 0, Ti < (d) = 0, > 







(e) r,. < 0, Ti < (f) > 0, Ti > (g) > 0, Ti < (h) < 0, Ti > 



Figure 1. All possible flowfields at the singular point 2 = associated with the dynamical 

system |lTJ. 

A collection of N of these point singularities, each located at z = zp{t), /3 = 
1, ■■■,N, by linear superposition, produces the field: 

1 ^ r 

i* = — zit) = x{t) + iyit)ec, r^ec. (1.6) 
0=1 ^ 

Then, if we advect each by the velocity field generated by all the other^ we arrive 
at the complex dynamical system: 

1 ^ r 

ia = i^J2' 1 ' z^{t) = Xc.{t) + iyo.{t)&^: r^eC, (1.7) 

ZTTi Zq, Zjj 

where ' indicates that (3 ^ a. In this paper we characterize all fixed equilibria 
of (1.7), namely solutions for which i*(i) = 0. For this, we have the N coupled 
equations: 

Y,' ^=0, (« = l,...iV), (1.8) 



f One might characterize this dynamical assumption by saying that each singularity 'goes with 
the flow'. 
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where we are interested in positions and strengths Ta for which this nonUn- 



ear algebraic system is satisfied. Since Eqn (1.8) is linear in the F's, it can more 
productively be written in matrix form 

Af = (1.9) 

where A E (£Nxn evidently a skew-symmetric matrix A — ~A^, with entries 
[o-aa] — 0, [cta^] — ^ — —[apa]- Wc Call A the configuration matrix associated 



with the interacting particle system ( 1.7 1. The collection of points {zi(0), Z2{0), zn(0)} 



in the complex plane is called the configuration. From (1.9 1, we can conclude that 
the points Za are in a fixed equilibrium configuration if det(A) = 0, i.e. there is at 
least one zero eigenvalue of A. If the corresponding eigenvector is real, the configu- 
ration is made up of point-vortices. If it is imaginary, it is made up of sources and 
sinks. If it is complex, it is made up of spiral sources and sinks. Notice also that if 
= 0, then one can prove that -gpf^ — for any n. It follows that: 

Theorem 1. For a given configuration of N points {zi, Z2, zn} in the complex 
plane, there exists a set of singularity strengths T for which the configuration is 



a fixed equilibrium solution of the dynamical system {1.7) iff A has a kernel, or 
equivalently, if there is at least one zero eigenvalue of A. If the nullspace dimension 
of A is one, i.e. there is only one zero eigenvalue, the choice of T is unique (up 
to a multiplicative constant). If the nullspace dimension is greater than one, the 
choice of T is not unique and can be any linear combination of the basis elements 
of null(A). 

The equilibria we consider in this paper all have one-dimensional nuUspaces and 
odd N . The more delicate cases of equilibria with higher dimensional nuUspaces 
and even N are deferred to a separate study. We mention here work of Campbell 
& Kadtke (1987) and Kadtke & Campbell (1987) in which a different technique is 



described to find stationary solutions to (1.7) 



2. General properties of the configuration matrix 

Since A is skew-symmetric, it follows that 

det(A) = det(-A^) = (-1)^ det(A^) ^ Aci[A^). (2.1) 

Hence, for N odd, we have — det(^^) — dei{A^), which implies det(^-^) = 0. 

Theorem 2. When N is odd, A always has at least one zero eigenvalue, hence 
for any configuration there exists a choice F G C for which the system is a fixed 
equilibrium. 

When N is even, there may or may not be a fixed equilibrium, depending on whether 
or not A has a non-trivial nullspace. In general, we would like to determine a basis 
set for the nullspace of A for a given configuration, i.e. the set of all strengths for 
which a given configuration remains fixed. Other important general properties of 
skew-symmetric matrices are listed below: 

1. The eigenvalues always come in pairs ±A. If N is odd, there is one unpaired 
eigenvalue that is zero. 
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2. If N is even, det(^) = PfiA)^ > 0, where Pf is the Pfaffian. 

3. Real skew-symmetric matrices have pure imaginary eigenvalues. 

Recall that every matrix can be written as the sum of a Hcrmitian matrix 
{B = B^) and a skew-Hermitian matrix (C = — C^). To see this, notice 

^^i(A + At) + i(A-At). (2.2) 

Here, B = ^{A + A^) = B^ and C = \{A- A^) = -C^ A matrix is normal 
if AA^ = A, otherwise it is non-normal. If we calculate AA^ — A^ A^ where 
A = B + C as, above, then it is easy to see that 

AA"^ - A'^A = 2(CB-BC). (2.3) 

Therefore, if _B = or C = 0, A is normal. 

Theorem 3. All Hermitian or skew-Hermitian matrices are normal. 



The generic configuration matrix A arising from (1.91 is, however, non-normal 



(a) Spectral decomposition of normal and non-normal matrices 
For normal matrices, the following spectral-decomposition holds: 
Theorem A. A is a normal matrix A is unitarily diagonalizable, i.e. 

A = QAQt (2.4) 

where Q is unitary. 

Here, the columns of Q are the N linearly independent eigenvectors of A that 
can be made mutually orthogonal. The matrix A is a diagonal matrix with the N 
eigenvalues down the diagonal. See Golub and Van Loan (1996) for details. 



In general, however, for the system of interacting particles governed by (1.8) 
(1.9 1, A G C^^^ will be a non-normal matrix. The most comprehensive decompo 



sition of A in this case is the singular value decomposition (Golub and Van Loan 
(1996), Trefethen and Bau (1997)). It is a factorization that greatly generalizes the 
spectral decomposition of a normal matrix, and it is available for any matrix. 

The TV singular values, a^"^^ (i = 1, . . . N), of A, are non-negative real numbers 
that satisfy 

Av«=a(*'u«; A^u«=a«v(^), (2.5) 

where u^'^ e and v*^*^ € C^. The vector u^*' is called the left-singular vector 
associated with a^-^\ while v^*' is the right-singular vector. In terms of these, the 
matrix A has the factorization 



A = UT.V^ = f^^*^u«v(*)^, (fc < N) (2.6) 

i=l 

where U e C^^^ in unitary, V G C^^^ is unitary, and S G M^^^ is diagonal. 



(2.6 1 is the non-normal analogue of the spectral decomposition formula (2.4) where 
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the summation term on the right hand side gives an (optimal) representation of A as 
a hnear combination of rank-one matrices with weightings governed by the singular 
values ordered from largest to smallest. Here, the rank of A is A:. The columns of 
U are the left-singular vectors u^*\ while the columns of V are the right-singular 
vectors v^^'^. The matrix E is given by: 



^(1) ... 



oNxN 



(2.7) 



AN) 



The singular values can be ordered so that cr*-^-' > ct*-^'' > . . . > cr^^-* > and one 



or more may be zero. As is evident from multiplying the first equation in (2.5) by 
and the second by A, 



(AU - aW2)^W ^ Q. (^^t _ cr«2)^W ^ 



(2.8) 



the singular values squared are the eigenvalues of the covariance matrices A^ A or 
AA\ which have the same eigenvalue structure, while the left-singular vectors u^^'^ 
are the eigenv ecto rs of AA\ and the right-singular vectors v*^') are the eigenvectors 
of At A. From (2.5), we also note that the right singular v ecto rs v^*' corresponding to 
o-(*) = form a basis for the nuUspace of A. Because of ( 1.9 ), we seek configuration 
matrices with one or more singular values that are zero. 



3. Collinear equilibria 

For the special case in which all the particles lie on a straight line, there is no loss 
in assuming — E M. Then A € M.^^^ , A is a normal skew-symmetric matrix, 
and the eigenvalues are pure imaginary. As an example, consider the collinear case 
N = 3. Let the particle positions be xi < X2 < x^, with corresponding strengths 
Ti, Fs. The A matrix is then given by 



A = 



X2—X1 
1 



3:1-2:3 
1 



(3.1) 



Since N is odd, we have Act{A) = 0. The other two eigenvalues are given by: 



\123 



±i, 



{x2-xi)^ (2:3 -2:2)2 {X3-Xiy 



(3.2) 



which is invariant under cyclic permutations of the indices (A123 = A231 = A312). 
We can scale the length of the configuration so that the distance between xi and 
0:3 is one, hence without loss of generality, let xi = 0,X2 = x,X3 = 1. The other 
two eigenvalues are then given by the formula: 



A — ±i\ 



'(1 



2:2(1 - xy 



(3.3) 
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It is easy to see that the numerator has no roots in the interval (0, 1), hence the 
nuUspace dimension of A is one. The nuUspace vector is then given (uniquely up to 
multiplicative constant) by: 



r = 



X3-X1 

X'j — Xi 



(3.4) 



For the special symmetric case X3 — xi = 1, — X2 ~ 1/2, X2 ~ xi — 1/2, we have 
Fi = 1, F2 = —1/2, Fs = 1. We show this case in figure [2] along with the separatrices 
associated with the corresponding flowfield generated by the singularities. Since the 
sum of the strengths of the three vortices is Fi+F2 + F3 = l — 1/2 + 1 = 3/2, the 
far field is that of a point vortex of strength F = 3/2. Interestingly, for the coUinear 
cases, since A is real, the nuUspace vector is either real, or if multiplied by i, is 
pure imaginary. Hence, each collinear configuration of point vortices obtained with 
a given F S M is also a collinear configuration of sources/sinks with corresponding 
strengths given by iF. The corresponding streamline pattern for the source/sink 
configuration, as shown in the dashed curves of figure [2] is the orthogonal comple- 
ment of the curves corresponding to the point vortex case. 




Figure 2. N = 3 evenly distributed point vortices on a line with strengths 
Fi = l,r2 = — |,r3 = 1, in equilibrium. The far field is that of a point vortex at the cen- 
ter-of-vorticity of the system. Solid streamline pattern is for point vortices, dashed streamline 
pattern is for source/sink system. The patterns are orthogonal. 



For N even, we cannot say a priori whether or not det{A) 
N — 2 shows. For this, the A matrix is 



as the case for 



L X2—X1 



X1—X2 







[ 







(3.5) 
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Figure 3. N = 7 evenly distributed point vortices on a line. The far field is that of a point vortex at 
the center-of-vorticity of the system. Because of the symmetry of the spacing, the vortex strengths 
are symmetric about the central point which also corresponds to the center-of-vorticity. 




The eigenvalues are A = ±i/d, hence there is no equihbrium (except in the hmit 
d ~^ oo). 

We show in figures [3] and |4] two representative examples of coUinear fixed point 
vortex equilibria for N — 7, along with their corresponding global streamline pat- 
terns. In figure |3] we deposit seven evenly spaced points on a line and solve for the 
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nuUspace vector to obtain the singularity strengths (ordered from left to right) 

f = (1.0000,-0.5536,0.9212,-0.5797,0.9212,-0.5536,1.0000), (3.6) 
= 2.1555. (3.7) 

a 

Because of the even spacing, the strengths are symmetric about the central point 
(Fi = Fy, r2 = Fg, Fa = F5), which is also the location of the center-of-vorticity 
Y!il^=i ^aXa- Figure [4] shows a fixed equilibrium corresponding to seven points ran- 
domly placed on a line. The nuUspace vector for this case is (ordered from left to 
right) 

f = (1.0000,-0.5071,0.5342,-0.4007,0.2815,-0.2505,1.0743), (3.8) 
^F„ = 1.7317. (3.9) 

a 

In both cases, the singularities are all point vortices (or source/sink systems) hence 
are examples of collinear equilibria such as those discussed in Aref (2007a, 2007b, 
2009) and Aref et. al. (2003) where the strengths are typically chosen as equal. 
The streamline pattern at infinity in both cases is that of a single point vortex of 
strength X]I=i Fq 7^ located at the center of vorticity ^a^a- 



4. Triangular equilibria 

The case iV = 3 is somewhat special and worth treating separately. Given any three 
points {zi, 22,23} in the complex plane, the corresponding configuration matrix A 
is: 



A 



21—2:2 21—23 
n 1 



22-21 
1 



22-23 




(4.1) 



23—21 23—22 

There is no loss of generality in choosing two of the points along the real axis, one 
at the origin of our coordinate system, the other at x = I. Hence we set zi = 0, 
Z2 — 1, and we let 23 = z. Then A is written much more simply: 



A = 



-1 

1 



1 

2-1 



1 

1-2 




(4.2) 



Since N is odd, one of the eigenvalues of A is zero. The other two are given by: 



(i-zy 



1 = ±i^ 



'(1-2 + 22)2 
22(1-2)2 



(4.3) 



When the numerator is not zero, the nuUspace dimension is one and it is easy to 
see that the nuUspace of A is given by: 



F = 



1 

2-1 
1 



(4.4) 
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However, the numerator is zero at the points: 

TTi Stti 

^ = exp(y),exp( — ), (4.5) 

at which 3?z = 3z ~ This forms an equilateral triangle in which case the 

nullspace dimension is three. We have thus proven the following: 

Theorem 5. For three point vortices, or for three sources/sinks, the only fixed 
equilibria are collinear. In this case, the nullspace dimension of A is one and is 



given by {4-4^ ■ For the equilateral triangle configuration, the nullspace dimension is 
three. 

We show a fixed equilibrium equilateral triangle state in figure [5] along with 
the corresponding streamline pattern. Figures [6j [7j [8] show examples of = 3 
triangular states that are not equilateral. 




Figure 5. = 3 equilateral triangle configuration with corresponding streamline pattern. The 
strengths are given by Ti = 1.0000, r2 = -0.5000 + 0.8660i, Ta = -0.5000 + 0.8660i. 



5. Equilibria along prescribed curves 

We now ask a more general and interesting question. Given any curve in the complex 
plane, if we distribute points {z^}, (a = 1, N) along the curve, is it possible to 
find a strength vector F so that the configuration is fixed? The answer is yes, if N 
is odd, and sometimes, if A'" is even. 

Figures [9 - [15] show a collection of fixed equilibria along curves that we prescribe. 
First, figure 9] shows 7 points places randomly in the plane, with the singularity 
strengths obtained from the nullspace of A so that the system is an equilibrium. The 
strengths are given by: f = (1.0000, -0.7958 + 1.0089i, -1.3563 - 0.4012i, 0.0297 + 
0.1594i, 0.9155 + 0.3458Z, -2.0504-0.8776i, -0.1935- 1.0802i)^ with the sum given 
by —2.4508 — 0.8449i. Thus, the far field is that of a spiral-sink configuration. Figure 
[lOlshows the case of = 7 points distributed evenly around a circle. The nullspace 
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Figure 6. N = 3 non-equilateral triangular state with corresponding streamline pattern. The 
strengths are given by Ti = 1.0000, r2 = 0.3077 + 0.4615i, Ts = -0.3200 - 0.2400j. 




Figure 7. N = 3 isosceles triangle state. The strengths are given by Fi = 1.0000, 
r2 = -1.6000 + O.SOOOi, Ts = -1.6000 - O.SOOOi. 

vector is given by f = (1.0000,-0.9010 + 0.4339^,0.6235 - 0.7818^,-0.2225 + 
0.9749i, -0.2225-0.9749i,0.6235 + 0.7818i, -0.9010-0.4339^)^. For this very sym- 
metric case, the sum of the strengths is zero, hence in a sense, the far field vanishes. 
Figure [TT] shows the case of N = 7 points placed at random positions on a cir- 
cle. Here, the nullspace vector is given by f = (1.0000, -0.6342 + 0.4086i, 0.3699 - 
0.5929i, -0.1501 0.6135i, -0.2483 - 0.9884i, 0.2901 -I- 0.3056i, -0.3595 - 0.2686ii)^ 
The random placement of points breaks the symmetry of the previous case and 
the sum of strengths is given by 0.2649 — 0.5222i which corresponds to a spiral- 
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sink. In figures [12] and [13] wc show equilibrium distribution of points along a 
curve we call a 'flower-petal', given by the formula r{9) = cos(20), < 9 < 27r. 
In figure [12] we distribute them evenly on the curve, while in figure 13 we dis- 
tribute them randomly. The particle strengths from the configuration in figure [12] 
are f = (1.0000, 0.1824 + 0.1498i, -0.9892 - 0.9103i, -0.1378 - 0.5333i, -0.1378 + 
0.5333^,-0.9892 + 0.9103^,0.1824 - 0.1498i)^ with sum equaling -0.8892 corre- 
sponding to a far field point vortex. Figure [13] shows particles distributed randomly 
on the same flower-petal curve. Here, the particle strengths are F = (1.0000, 0.2094— 
0.4071?, -0.3009-H0.3003i, 0.0404-0. 2864i, -0.1779-h0.2773z, 0.4236-h0.8052i, -0.4702- 
0.3304i)^, with sum given by .7244 -f .3589z. Hence the far field corresponds to a 
source-spiral. 

The last two configurations, shown in figures [14] and [15] are equilibria distributed 
along figure eight curves, given by the formulas r{9) — cos'^{9), < 9 < 2tt. 
In figure [14] we distribute the points evenly around the curve, which gives rise to 
strengths F = (1.0000, -0.2734-h0.5350i, 0.0239-0.2080i, 0.1063-0.0517i, 0.1063-h 
0.0517i, 0.0239 + 0.2080i, -0.2734 - 0.5350i)^, whose sum is .7136, thus a far field 
point vortex. In contrast , w hen the points are distributed randomly around the 
same curve, as in figure 15 the strengths are given by F = (1.0000,-0.1054 -f 
0.5724i, -0.0174-0.4587i, 0.9208+1. 2450i, -0.0460-0.4577i, -0.5292-1-0. 237H, -0.2543- 
0.0921i)^, with sum equaling .9685 -I- 1.0460i, hence a far field source-spiral. 



6. Classification of equilibria in terms of the singular 

spectrum 

Here we describe how to use the non-zero singular spectrum of A to classify the 
equilibrium when A has a kernel. Let cr'-'^ i = 1,...,A: < N denote the non-zero 
singular values of the configuration matrix A, arranged in descending order cr^^) > 
c*-^^ > ... > a'^^'> > 0. First we normalize each of the singular values so that they 
sum to one: 
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I I — _ I I I I 

Figure 9. Fixed equilibrium for seven points placed at random locations in the plane. The far 
field is a spiral-sink (figure 1(e)) with since J^^a = —2.4508 — 0.8449i. 



k 

^(«) (6.1) 

Then 

k 

^a« = l, (6.2) 

i=l 

and the string of k numbers arranged from largest to smallest: {a^^\a''^\ ...,a'^'^^) 
is the 'spectral representation' of the equilibrium. The rate at which they decay 
from largest to smallest is encoded in a scalar quantity called the Shannon entropy, 
S, of the matrix (see Shannon (1948) and more recent discussions associated with 
vortex lattices in Newton and Chamoun (2009)): 



S 



^aWlo, 



'go- 



(i) 



(6.3) 



With this representation, spectra that drop off rapidly from highest to lowest, are 
'low-entropy equilibria', whereas those that drop off slowly (even distribution of 
normalized singular values) are 'high-entropy equilibria'. Note that from the rep- 



resentation (2.6), low-entropy equilibria have configuration matrix representations 
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Figure 10. N = 7 evenly distributed points on a circle (dashed curve) in equilibrium. Because of 
the symmetry of the configuration, ^ Ta = 0, hence the far-field vanishes. 



that are dominated in size by a small number of terms, whereas the configuration 
matrices of high-entropy equilibria equilibria have terms that are more equal in size. 
See Newton and Chamoun (2009) for more detailed discussions in the context of 
relative equilibrium configurations, and the original report of Shannon (1948) which 
has illuminating discussions of entropy, information content, and its interpretations 
with respect to randomness. As an example of the normalized spectral distribution 
associated with the figure eight equilibrium shown in figure |15[ we show in figure 



16 the 7 singular values (including the zero one). The fact that they are grouped 



in pairs follows from the skew-symmetry of A which implies that the eigenvalues 
come in pairs ±A. Since the singular values are the squares of the eigenvalues, it 
follows that there are two of each of the non-zero ones. 

Tables [T]-[5]show the complete singular spectrum for all the equilibria considered 
in this paper. A common measure of 'robustness' associated with the configuration 
matrix, hence the equilibrium, is the size of the 'spectral gap' as measured by 
the size of the smallest non-zero singular value. From Table [2] the coUinear state 
with points distributed randomly and the figure-eight state with points distributed 
evenly (Table |4]) are the least robust in that their smallest non-zero singular values 
are closest to zero. 
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Figure 11. N = 7 randomly distributed particles on a circle (dashed curve) in equilibrium along 
with the corresponding streamline pattern. The far field streamline pattern is that of a spiral-sink 
(figure 1(g)) since J2^a = 0.2649 - 0.5222i. 



7. Discussion 



In this paper we describe a new method for finding and classifying fixed equilibrium 
distributions of point singularities of source/sink, point vortex, or spiral source/sink 
type in the complex plane under the dynamical assumption that each point 'goes 
with the fiow'. This includes configurations placed at random points in the plane, at 
prescribed points, or lying along prescribed curves. This last situation is reminiscent 
of a classical technique for enforcing boundary conditions along arbitrarily shaped 
boundaries embedded in fiuid flows. These techniques are generically referred to 
as singularity distribution methods. See, for example, Katz and Plotkin (2001) 
and Cortez (1996, 2000) for applications and discussions of these methods in the 
context of potential flow, hence inviscid boundary conditions, and Cortez (2001) in 
the context of Stokes flow, hence viscous boundary conditions. For these problems, 
there is generally no associated evolution equation for the interacting singularities 
which discretize the boundary, as in (1.7) for us. Their positions are fixed to lie 



along the given boundary, and the strengths are then judiciously chosen to enforce 
the relevant inviscid or viscous boundary conditions. As in Cortez (2001), it would 
be of interest to 'regularize' the point singularities (1.1) and ask if the methods in 
this paper can be extended to smoothed out singularities, as would a more complete 
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Figure 12. N = 7 evenly distributed particles in equilibrium on the curve r{6) = cos(20) (dashed 
curve) along with the corresponding streamline pattern. The far field corresponds to a point vortex 
since J^Td = -0.8892. 

analysis of the 'pscudo-spcctrum' associated with the configuration matrices A, as 
discussed Trefethen and Bau (1997). 
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Figure 13. N = 7 randomly distributed particles in equilibrium on the curve 
r{9) = cos(2S) (dashed curve). The far field corresponds to a source-spiral (figure 1(f)) since 
= 0.7244 + 0.35891. 
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Figure 16. Singular spectrum (normalized) for N = 7 particles placed randomly along a figure 
eight planar curve (i.e. the equilibrium configuration shown in figure [Tsl l . Singular values are 
grouped in pairs (except for the zero value) due to the skewr-symmetry of the configuration matrix. 



Configuration 


a (unormalized) 


a (normalized) 


Shannon entropy 




1.0000 


0.5000 


0.6931 


Equilateral 


1.0000 


0.5000 






0.00 


0.00 






1.0598 


0.5000 


0.6931 


Isosceles (acute) 


1.0598 


0.5000 






0.00 


0.00 






2.7203 


0.5000 


0.6931 


Isosceles (obtuse) 


2.7203 


0.5000 






0.00 


0.00 






1.2115 


0.5000 


0.6931 


Arbitrary triangle 


1.2115 


0.5000 






0.00 


0.00 





Table 1. Singular spectrum of triangular states (N = 3) 
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Configuration 


a (unormafized) 


a (normalized) 


Shannon entropy 




4.5000 


0.5000 


0.6931 


N = 3 


4.5000 


0.5000 






0.00 


0.00 






2.5249 


0.3214 


1.5237 


N = 7 (even) 


2.5249 


0.3214 






1.6831 


0.1428 






1.6831 


0.1428 






0.8420 


0.0357 






0.00 


0.00 






6.3408 


0.4457 


1.0723 


N = 7 (random) 


6.3408 


0.4457 






2.0969 


0.0487 






2.0969 


0.0487 






0.7062 


0.0055 






0.7062 


0.0055 






0.0000 


0.0000 





Table 2. Singular spectrum of collinear states (N = 3, 7) 



Configuration 


a (unormalizcd) 


a (normalized) 


Shannon entropy 




3.0000 


0.3214 


1.5236 


N = 7 (even) 


3.0000 


0.3214 






2.0000 


0.1429 






2.0000 


0.1429 






1.0000 


0.0357 






1.0000 


0.0357 






0.0000 


0.0000 






3.7954 


0.3363 


1.4700 


N = 7 (random) 


3.7954 


0.3363 






2.4250 


0.1373 






2.4250 


0.1373 






1.0631 


0.0264 






1.0631 


0.0264 






0.0000 


0.0000 





Table 3. Singular spectrum of circular states (N = 7) 
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Configuration 


(T (unormalizod) 


(T (normalized) 


Shannon entropy 




11.9630 


0.4664 


0.9651 


N = 7 (even) 


11.9630 


0.4664 






3.0001 


0.0293 






3.0001 


0.0293 






1.1454 


0.0043 






1.1454 


0.0043 






0.0000 


0.0000 






6.9337 


0.3465 


1.3929 


N = 7 (random) 


6.9337 


0.3465 






4.4357 


0.1418 






4.4357 


0.1418 






1.2769 


0.0117 






1.2769 


0.0117 






0.0000 


0.0000 





Table 4. Singular spectrum of figure eight states (N = 7) 



Configuration 


a (unormalized) 


a (normalized) 


Shannon entropy 




5.9438 


0.4447 


1.1034 


TV = 7 (even) 


5.9438 


0.4447 






1.8115 


0.0413 






1.8115 


0.0413 






1.0538 


0.0140 






1.0538 


0.0140 






0.0000 


0.0000 






8.0780 


0.3875 


1.3393 


N = 7 (random) 


8.0780 


0.3875 






3.8900 


0.0899 






3.8900 


0.0899 






1.9523 


0.0226 






1.9523 


0.0226 






0.0000 


0.0000 





Table 5. Singular spectrum of flower states (N = 7) 
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